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ABSTRACT 


A  prediction  method  for  estimating  wheel  sinkage,  landing 
gear  drag  and  aircraft  take-off  performance  on  clay  and  sand 
airfields,  of  various  bearing  strengths,  is  described.  Estimates 
of  take-off  distance,  wheel  drag  to  vertical  force  ratios  for  the 
full  ground  velocity  range,  and  the  wheel  sinkage  into  the  soils 
have  been  prepared  for  the  C-5  and  C-141  aircraft.  These  charts 
have  been  prepared  for  various  aircraft  operating  gross  weights 
and  selected  tire  pressures.  Results  indicate  that  the  C-5  airplane 
can  operate  successfully  at  a  gross  weight  of  571,000  pounds,  with 
95  PSI  tire  pressure,  from  soil  airfields  having  a  bearing  strength 
of  CBR  »  5.  The  take-off  distance  is  increased  by  205!  over  paved 
surfaces  and  the  maximum  wheel  sinkage  into  the  soil  is  about 
1.2  inches.  The  C-141  airplane  can  operate  at  a  gross  weight  of 
280,000  pounds,  with  100  PSI  tire  pressure,  from  soil  airfields 
having  a  bearing  strength  of  CBR  -  5.5.  The  increase  in  take-off 
distance  is  about  205!  greater  than  for  paved  surfaces  and  the 
maximum  wheel  sinkage  into  the  soil  is  about  1.6  inches. 

Three  FORTRAN  IV  computer  programs,  SOLDG,  S0LDG2  and  TAKOFF 
for  computing  landing  gear  drag,  have  been  prepared  and  are  listed 
in  the  report.  SOLDG  is  employed  when  only  the  weight  and  tire 
configuration  are  known  and  S0LDG2  require^  knowledge  of  the 
airplane  aerodynamic  lift  and  landing  gear  geometry.  Both  of  these 
programs  assume  constant  taxi  velocities.  TAKOFF  is  a  more 
sophisticated  program  which  requires  knowledge  of  the  airplane 
aerodynamic  lift,  landing  gear  geometry,  and  power  plant  thrust 
and  as  used  to  compute  take-off  distances.  Example  problems  using 
each  of  the  computer  programs  are  illustrated. 
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SECTION  I 


INTRODUCTION 


Operating  of  aircraft  from  minimum  length  runways  with  substandard 
strengths  has  been  given  increased  emphasis  in  view  of  present  opera¬ 
tional  requirement  estimates.  In  anticipation  of  forthcoming  requests 
for  methods  to  predict  loads  and  performance  of  aircraft  operating  in 
remote  areas  where  good,  long,  hard  surfaced  runways  would  not  be 
available,  a  program  was  initiated  in  1964  with  The  Boeing  Company  to 
provide  some  analytical  means  by  which  these  performance  estimates 
could  be  made.  This  study  has  recently  been  completed  and  as  such, 
provides  the  first  analytical  model  available  for  the  prediction  of 
aircraft  soil  drag,  sinkage  and  takeoff  performance  on  clay  and  sandy 
type  soils. 

Various  techniques  were  used  to  predict  soil  static  and  dynamic 
effects,  however,  the  method  which  was  finally  considered  to  give  the 
best  correlation  with  available  data  was  basically  one  developed  by  the 
U.S.  Army  Waterways  Experiment  Station  (WES)  modified  for  high  speed 
effects.  This  method  uses  empirical  sinkage  data  obtained  by  towing 
carts  with  various  loads  and  various  tire  sizes  over  especially 
prepared  soil  beds.  The  results  were  collapsed  i^to  a  dimensionless 
number  called  the  Mobility  Number  which  is  a  function  of  the  wheel, 
the  load  and  the  soil  str;  ngth  characteristics . 

The  basic  Mobility  Numbers  were  modified  by  Boeing  to  account  for 
dynamic  effects  not  covered  by  the  original  relationships.  They  were 
then  used  to  determine  analytically  the  performance  of  the  367-80 
aircraft  (KC-135  prototype)  and  compare  these  results  to  a  high- 
flotation  test  conducted  on  Harper's  dry  lake. 

The  purpose  of  this  technical  report  is  to  present  three  programs 
for  automating  the  methods  generated  by  Boeing  and  WES  and  to  calculate 
the  performance  capability  of  two  other  first  lire  cargo  type  aircraft 
(C-141  and  the  C-5)  while  operating  on  runways  of  varying  soil  strength. 

The  three  programs  cover  varying  levels  of  detail  which  may  be 
required.  The  first  is  used  to  calculate  drag  and  sinkage  for  wheels 
attached  to  one  strut  while  taxiing  at  constant  speed.  The  second  con¬ 
siders  the  complete  aircraft  taxiing  at  constant  velocity  and  the  last 
calculates  takeoff  performance  of  a  complete  aircraft  on  hard  surface 
(paved  runways)  and  on  clay  and  sandy  soils.  The  computer  programs  were 
tested  for  accuracy  by  analyzing  conditions  which  were  identical  to  those 
uaeu  by  Boeing  and  comparing  the  resuits  with  Boeing's  full  scale  air¬ 
plane  tests  at  Harper  Laxe.  These  results  were  identified  as  "FDDS 
Correlation." 

Finally,  a  discussion  is  given  of  the  limitations  of  the  method 
and  of  ways  by  which  improvements  may  be  achieved. 
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SECTION  II 


DESCRIPTION  OF  MODEL 


1.  Discussion  -  Although  Reference  1  gives  a  detailed  development  of 
the  analytical  wheel-soil  model,  a  summarized  description  is  given  here 
for  completeness. 

brag  loads  resulting  from  a  tire  rolling  on  a  soft,  yielding  soil 
can  be  considered  to  be  made  up  of  two  interrelated  components,  a  static 
(low  speed)  and  dynamic  (high  speed)  effect.  The  equations  used  in  the 
TAKGFF  program  will  be  derived  frou  basic  relationships.  The  other  two 
programs.  SGLDG2  and  SOLDG,  are  merely  formed  as  special  cases  of  the 
general  development. 

An  aircraft  accelerating  under  the  thrust  of  engines  obeys 
the  familiar  relation: 

IF  =  sa 


where 


IF  =  The  sum  of  all  external  forces  in  the  plane 
cf  the  motion 

b  =  Hass  of  the  body 

a  =  The  body’s  acceleration  resulting  from  the 
applied  forces. 

The  forces  making  up  the  IF  can  be  seen  in  Figure  1.  The  problem 
in  the  present  study  is  to  determine  the  drag  forces  associated  with 
the  landing  gear  (F  and  in  the  figure).  The  c*her  forces  will 

"N  "x 

be  discussed  later. 


Forces  Acting  on  an  Aircraft  While  Taxiing 
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FIGURE 


On  r avec  surfaces  these  forces  are  merely  the  coefficient  of  rolling 
friction  u  (usually  taken  as  0.02  or  0.03  for  paved  surfaces)  times  the 
normal  •:  rce  or  ground  reactions  R  and  R^.  On  soil,  however,  a  more 

difficult  problem  exists.  Figure  5  sl.-'ws  a  tire  traversing  a  soil 
surface.  As  one  «ould  expect,  the  sinkage  and  consequently  tha  force 
required  to  move  the  wheel  is  greatly  affected  by  the  speed  of  tee 
vehicle . 


FV 


FIGURE  2  Tire  Traversing  a  Soil  Surface 


As  mentioned  above,  two  drag  forces  make  up  the  total  drag  of  the  wheels. 

In  Figure  2,  D  is  the  low  speed  rolling  resistance  force  and  D  is  the 
K  I 

high  speed  impingement  drag  force.  These  forces  will  be  discussed 
separately,  although  some  interaction  exists  because  both  are  functions 
of  some  of  the  same  parameters. 

1.1  Soil  Strength  -  Before  discussing  the  details  of  these  drag 
forces,  a  few  words  are  in  order  to  describe  the  basis  for  determining 
soil  strength.  In  this  study  and  in  much  of  the  work  being  done  in 
soil  mechanics,  strength  is  measured  by  the  California  Bearing  Ratio 
(CBR)  which  is  a  measure  of  the  resistance  to  penetration  that  tne  soil 
exhibits  to  a  cylinder  of  three  sq.  in.  area.  That  is,  a  piston  or 
bearing  plate  is  subjected  to  an  increasing  load  and  deformations  are 
noted.  The  ratio  of  the  bearing  load  per  unit  deformation  to  the 
deformation  of  a  standard  soil  (California  limestone  gravel)  in  percent 
is  the  CBR.  Details  of  CBR  measurements  are  in  accordance  with  MIL - 
STi’-621.  Because  of  the  difficulty  in  making  CBR  measurements,  WES 
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developed  the  airfield  penetrometer  which  is  a  device  consisting  of  a 
30— degree  cone  with  base  area  of  0.196  square  inches  mounted  on 
the  bottom  of  a  graduated  rod.  When  the  penetrometer  Is  forced  into 
the  ground,  a  dial  gage  mounted  on  top  of  the  rod  gives  a  measure  of 
the  force  required  to  penetrate  the  soil.  The  reduced  readings  of  this 
device  are  known  as  Airfield  Index  (AI).  The  Cone  Index  (Cl)  which  is 
used  in  the  following  forme lations  is  based  on  an  identical  procedure 
as  that  of  AI  except  that  the  scale  and  cone  area  are  different.  The 
scale  used  in  Ci  measurement  (0  -  300)  represents  a  force  of  0  -  150 
pounds  which  results  from  a  cone  having  an  area  of  0.50  square  inches. 
Obtaining  CI  directly  requires  too  much  effort  for  hand  penetration 
thus,  the  Airfield  Index  requiring  about  one^fifth  the  effort  was 
adopted.  Conversion  from  AI  to  CI  is  merely  a  scale  change  and  is  a 
1  to  50  rrtio.  Conversion  from  CBR  to  AI  is  approximately  a  lto  1 
ratio,  however,  this  may  vary  widely  and  caution  should  be  used  in 
manipulating  CBR  and  CI  to  avoid  errors  in  conversion.  Figure  3 
shows  typical  plots  of  AI  and  CBR.  Notice  that  values  range  froir 
•  45tol,  to  1.4  to  1,  for  CBR  to  AT.  ratios.  It  would  be  advisable  to 
measure  CI  directly  and  avoid  the  conversion  altogether .  Boeing 
indicates  that  if  the  AI  and  CBR  are  known,  the  conversion  to  CI  can 
be  made  using  the  50 to  1  ratio  mentioned  above,  but  if  only  CBR  is 
known,  an  approximate  aud  widely  accepted  conversion  is  lto  40.  This 
value  will,  however,  underpredict  sinkages  as  given  by  the  present 
method,  and  a  more  realistic  ratio  of  about  lto  25  CBR  to  CI  should  be 
used.  In  the  present  studies,  an  assumption  is  made  that  CI  is  known, 
or  that  when  a  correlation  is  made  for  the  367-80  data  on  a  soil 
measured  in  CBR,  a  36  to  1  conversion  is  used  to  conform  with  the 
Boeing  assumptions. 


1.2  Sinkagf  -  One  other  parameter  which  must  be  discussed  before 
the  drag  terms  can  be  defined  in  detail  is  sinkage.  Sinkage  relation¬ 
ships  used  in  the  present  study  are  obtained  from  WES  empirical  data 
generated  by  towing  weighted  carts  through  soils  of  various  CI  using 
a  wide  variety  of  tire  geometry.  This  data  was  collapsed  Into  a  non- 
dimensional  term  called  "Mobility  Number". 

yl/2 

fi  -  1CI)B  D  {^t\  for  day  (1) 


and 


fej 


G(3  D) 


3/2 


ft) 


for  sand 


(2) 


where 


Rc,ft  *  Mobility  Numbers  for  clay  and  sand  respectively 
CI  ■  Cone  Index  (Psi) 


B  *  Width  of  Tire  (in.) 
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Site  A- 7 

Fort  Leonard  Wood,  MO. 


Site  U-8 

Fort  Campbell,  KY. 


Site  U-5 

Fort  Stewart,  GA. 


Site  U- 11 

Yuna  Proving  Ground,  Arizona 


Fort  Campbell,  KY. 


FIGURE  3  CBR  -  Airfield  Index  Correlation 


b 


D  =  Diameter  of  Tire  (in.) 

6  =  Tire  Deflection  Caused  by  on  a  rigid 

C  surface  (in.) 

=  Tire  Section  Height  (in.) 

Fy  =  Vertical  Force  on  a  Single  Tire  (lb.) 

G  =  Slope  of  Cl  vs  Depth  Curve  (Psi/in)  =  2  CI^^/B 

Using  these  equations ,  values  of  9  and  were  plotted  against  an 

incremental  sinkage  coefficient  Z/D  for  three  passes  over  the  same  path. 
Equations  were  then  generated  to  fit  these  data  for  a  range  of  mobility 
numbers.  The  resulting  curves  are  shown  on  Figures  4  and  5  for  the 
first  two  passes.  For  the  present  study,  only  first  nass  sinkage  data 
will  be  used. 


Z  _  0.1208 

D  *  S  -  0.9468 
c 


0.0095  for  clay 


(3) 


and 


Z  _  0.3439 
D  E  ~0.6239 

s 


0.0017  for  sand 


(4) 


These  equations  are  valid  for  fl  from  2.0  to  11.0  and  9  from  3.5  to 

c  s 

82.5, respectively.  Accuracy  is  affected  most  in  the  lower  range  as  can 
be  seen  from  the  figures. 

Since  these  data  were  obtained  at  a  very  slow  tewing  velocity,  some 
consideration  must  be  given  to  the  dynamic  effects  if  sinkages  at  higher 
velocities  are  to  be  valid.  This  was  done  by  fitting  an  equation 
through  experimental  data  to  determine  a  factor  which  when  multiplied 
times  the  slow  speed  sinkage  gives  th-:  high  speed  dynamic  sinkage.  The 
equation  for  this  dynamic  factor  is: 

1.0  +  1.34e"1,27tp  (5) 


Dynamic  Factor 

Pulse  time  or  time  of  application  of  deforming 
force  =  1  (sec) 

V~ 

Tire  Footprint  Length  (in.) 

Ground  Velocity  (in. /sec) 
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DF  = 

where 

DF  = 

t  = 
P 

where 

1  = 
t 

V  = 


0.14 


2nd  Pass 


Figuro  4  Incremental  Siukage  Coefficient  fur  Clay 


? 


INCREMENTAL  SINKAGE  COEFFICIENT,  Z/D 


i 


Figure  5  Incremental  Sinkage  Coefficients  for  Sand 
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10.0 


1.00 

PULSE  TIME,  t  -SECONDS 
P 


as  a  Function  of  Pulse  Time 
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In  the  equation  for  pulse  time,  the  footprint  is  a  function  of  the 
axle  displacement  and  becomes: 

lt  -  2  [  D.ZA  -  ZA2]  1/2  (6) 

where 


D  -  Tire  Diameter  (in.) 


ZA  -  Axle  Displacement  ■  6t  +  Z  or  tire  deflection 

The  variation  of  dynamic  factor  with  pulse  lime  is  shown  in  Figure  6. 
Using  the  dynamic  factor,  an  effective  soil  strength  is  determined  by 
modifying  the  Cl  or  G  of  the  soil  as  follows: 

Cl  -  DfjCI  (7) 

?ff  1.6 


G 


eff 


DFjG 

1.6 


(8) 


In  equations  7  and  8,  the  1.6  is  used  to  correct  the  basic  WES  data 
for  the  very  slow  towing  speeds.  This  addition  of  the  dynamic  effect 
now  allows  us  to  consider  the  mobility  number  as  a  dynamic  mobility 
number  and  gives  the  corrected  sinkage  coefficients  from  Figures  4  and  5. 
The  equation  becomes: 


and 


where 


Z'  _  0.1208 

D  “  0*  -  0.9468 
c 


0.0095  for  clay 


(9) 


V  0.3429 _ 

D  “  O'  -  0.6239 
s 


0.0017  for  sand 


(10) 


.1 '  -  Dynamic  Sinkage  (in.) 

D  -  Tire  Diameter  (in.) 

O'  »  Dynamic  Clay  Mobility  Number  - 

O'  *  Dynamic  Sand  Mobility  Number  ■  ^  .  0 

s  1 .6  s 


Calculation  of  dynamic  sinkage  is  an  iterative  process  since  the  dynamic 
factor  is  a  function  of  the  soil  sinkage.  The  present  program  assumes 
an  initial  sinkage  value  and  iterate*  continuously  updating  the  assumed 
value  until  the  calculated  sinkage  is  within  a  specified  percentage  of 
the  assumed  sinkage.  The  process  converges  quite  rapidly  to  the  actual 
s inkage . 
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Once  the  sinkage  values  are  determined,  the  drag  components  can  be 
determined.  The  next  section  discusses  the  procedure  for  determining 
these  parameters. 

2.  Soil  Drag 

In  the  previous  discussion,  two  types  of  soil  drag  were  briefly 
mentioned  -  the  low  speed  or  rolling  resistance  drag  and  the  impingement 
or  high  speed  drag.  This  section  will  discuss  each  type  individually 
although  an  interaction  can  be  observed. 


2.1  Low  Speed  Drag  -  The  drag  contributed  by  low  speed  sinkage, 
friction  or  rolling  resistance  considerations  can  be  compared  to  the 
classic  friction  drag  calculated  by  multiplying  the  normal  force  by 
some  friction  coefficients  (usually  0.02  or  0.03  for  paved  surfaces) 
as  seen  in  Equation  11. 


(ID 


The  coefficient  p  is  then  the  ratio  of  the  arag  force  over  the  vertical 
force.  The  soil  drag  was  determined  for  the  present  method  from  addi¬ 
tional  WES  empirical  drag  ratio  or  drag  coefficient  data  gathered  from 
low  speed  tow  tests  of  configurations  given  various  mobility  numbers. 


Curves  fitted  through  these  data  are  shown  as  Figures  7  and  8 
for  clay  and  sand.  These  curves  can  be  represented  by  the  following 
equations.  For  clay  soil  on  the  first  and  second  pass: 


D 

F- 

V 


0.4072 
C,  -  0.8713 


0.0206 


(12) 


and  for  sandy  soil  on  the  first  yaf,s: 


V 

and  on  the  second  pass : 


0.6490 

-  2.2222 


+  0.0366 


c 

“0 


* 


"2 


0.7085 
Gs  -  2.1739 


+  0.0195 


(13) 


(14) 


Dynamic  effects  are  again  considered  by  applying  the  dynamic  factor  as 
before  and  substituting  the  dynamic  mobility  number  for  w  and  0  ,  in 
equation*  12  and  13.  8 

2.2  High  Speed  Drag  -  The  high  speed  drag  conti i but ion  is  deter¬ 
mined  by  assuming  that  soil  drag  is  analogous  to  other  types  of  fluid 
flow  in  that  it  varies  proportionally  to  the  frontal  or  flat  plate  3rea 
which  the  tire  presents  to  the  soil.  This  drag  is  formulated  using  the 
familiar  equation 


D 


I 


CD^S 


(15) 


! 


DRAG  RATIO,  u* 


Figure  7 


Drag  Ratio  for  Light  Loads  on  Clay 


where 


D  =  Impingement  Drag  rlb.) 

I 

CD  =  Lxpinr-^ent  Drag  Coefficient 

q  =  Dynamic  Fressure  of  Soil  (Psi)  =  l/2p  V~ 

s 

=  Soil  Density  (lb/ sec*-/ m-*) 

V  =  Forward  Velocity  (ips) 

S  =  Tire  Frontal  Area  =  B*?' 

An  illustration  shoving  the  area  over  which  the  drag  is  applied  is 
given  in  Figure  9a.  Using  this  basic  assiaapticr. ,  a  series  >i  empirical 
C..  values  were  determined  such  that  the  equation  is  satisfied.  Tie 
UI 

drag  indicated  by  Equation  15  reaches  a  limit  at  some  velocity  bv  a 
lifting  effect  cauo^i  by  the  tire  planing  over  the  soil.  The  velocity 
vhere  planing  is  experienced  is  calculated  by  assuming  that  the  lift 
experienced  by  the  tire  follows  a  similar  celatit oship: 


where 


c-r  ~  C^qS 


(16) 


=  Lift  due  to  impingement  drag  (lb)  - 

Vertical  Tire  Force  F,. 

« 

=  Soil  Lift  Coefficient 


q  =  Dynamic  Scil  Pressure  (psi)  =  1/2"  V~ 
A.  =  Area  overvhica  F  is  distributed 

L  ■' 


2  »  4 

"■  s  =  Soil  Density  (lb/sec  /in  ) 

^  =  Tire  Planing  Velocity  (ips) 
Rea r r a ng f ag  Foliation  ±6  and  solving  for  V  ye  haver 


V 


P 


A  o  CT 
t  s  Lj. 


(17) 


14 


An  illustration  shoving  the  area  over  vhich  the  lift  force  is  assumed 
to  act  is  given  in  Figure  9b.  Experimental  values  of  Cj  to  cause 

planing  for  various  soil  strengths  were  determined  from  Boeing's 
Harper  Lake  test  series  and  are  shovn  in  Figure  10.  The  planing 
velocity  effect  on  the  impingement  drag  coefficient  is  shown  in 
Figure  11.  More  discussion  of  these  curves  is  given  in  Section  7. 

2.3  Total  Soil  Drag  -  The  total  soil  drag  can  new  be  determined 
from  the  components  given  by  paragraph  2.1  and  2.2  by  applying  them  to 
each  tire  of  the  vehicle  and  swing  up  to  obtain  the  total  drag.  The 
only  exception  to  this  is  when  the  geometry  is  such  that  a  tandem  set 
exists.  Then  the  drag  of  the  leading  tire  is  calculated  and  the  drag 
of  the  trailing  tire  is  assumed  to  be  a  percentage  of  the  leading  tire. 
This  slight  simplification  has  proved  to  be  a  reasonable  representation 
cf  the  drag  values.  The  total  drag  for  a  set  of  main  landing  gear 
then  becomes 


for  tandem  wheel  sets  and 


(18) 


(19; 


for  single  or  dual  wheel  sets. 

In  Equations  18  and  19: 

F  =  Total  Main  Landing  Gear  Drag  (lb) 

.  i 

Fj)  =  Soil  drag  for  a  single  set  of  Landes  Wheels 
T  (lb) 

N  =  Number  of  Tandem  Wheel  Sets  per  Side 
F„  =  Soil  Drag  for  a  Single  Tire  (lb: 


Sy=  Number  of 


The  terms  F  and  F  are  formed 
D  Dc 

T  S 


Tires  per  side. 

from  the  components  of  low  and  high  speed 


drag  mentioned  in  the  preceding  paragrenhs  and  are  given  by  equations 
20  and  21. 


15 


Experimental  Tire  Lift  Cocffi 


IMPINGEMENT  DRAG  COEFFICIENT 


Q 

o 


GROUND  VELOCITY  ,  V/Vp 
PLANING  VELOCITY 


Figure  11  Impingement  Drag  Coefficient 
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(20) 


r  *  (  *  +  *)  f„  +  Dr  (1  +  TWF)  for  sand 

\  vw i  »2  V 


F  -  2p*  Fv  +  Dj  (1  +  TWF)  for  clay 
°T 

jj*  pfc  “  Rolling  Resistance  Coefficient  from 

’  ’  Equations  L2,  13  and  14. 

F  *  Net  downward  vertical  load. 


D  =  Impingement  drag  from  Equation  15. 

TWF  ■  Tandem  wheel  factor  for  simplified  rear 

wheel  sinkage  and  drag,  (0.4  for  clay  soil 
and  1.0  for  sand  soil) 

Similarly,  for  a  single  or  multiple  side  by  side  nose  gear  wheel, 

F  =  (u*  Fv  +  DI)  NW  for  claM 
°S 
or 

F  «  (y*  Fy  +  Dj)  Nw  for  sand. 

DS 

3.  Formulation  of  Ground  Acceleration  Equation  -  Having  defined  the 
soil  drag  force,  we  can  formulate  the  ground  acceleration  equation  which 
when  integrated  will  yield  ground  distance.  Referring  to  Figure  1 
and  summing  forces  in  a  horizontal  direction,  we  have 

IF  =  ma 

or  solving  for  a  and  making  substitutions  for  £F 


a  =  8  T  -  D  -  (F  +  F  )  -  W- 
V  a  DN  °M 


where 


a  -  Aircraft  ground  acceleration  (ft/sec2) 
g  =  Acceleration  due  to  gravity  (ft/sec2) 

W  -  Gross  weight  (lb) 

T  *  Total  engine  thrust  (lb) 

*=  Aerodynamic  drag  (lb)  -  C^qA^ 
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Fp  ■  Soil  drag  from  nose  tires  (lb) 

N 

Fd  “  Soil  drag  from  main  tires  (lb) 

M 

<P  ■  Runway  slope  (Rad) 

Integration  of  equation  22  in  a  step-wise  manner  using  a  finite 
difference  technique  yields  a  time  history  of  parameters  such  as  soil 
drag,  velocity,  acceleration,  sinkage  and  takeoff  distance  cs  a 
function  of  time. 

Combining  these  results  with  equations  of  motion  for  the  landing 
gear  and  aircraft  structure  will  permit  studies  of  dynamic,  response 
due  to  traversing  soil  airfields  with  non-level  profiles.  Limited 
studies  of  this  type  have  shown  a  promising  step  toward  a  solution 
to  soft  field  dynamic  response  prediction  problems. 
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SECTION  III 


COMPUTER  PROGRAM  DESCRIPTION 


Tue  equations,  theory  and  assumptions  of  Section  II  have  been  used 
to  generate  three  computer  programs  which  may  be  used  to  calculate 
soil  sinkage,  drag  and  takeoff  performance  on  a  clay  or  sandy  soil. 

Three  programs  were  generated  so  that  the  problem  of  landing  gear  drag 
could  be  considered  when  varying  amounts  of  detailed  information  about 
the  vehicle  are  known. 

The  simplest  program  designated  SOLDQ  is  designed  to  study  soil 
drag  when  only  the  tire  configuration  and  weight  are  known.  The  program 
assumes  sufficient  power  to  maintain  a  constant  taxi  velocity  and  zero 
aerodynamic  lift  and  drag.  The  second  program  designated  S0LDG2,  also 
assumes  a  constant  taxi  velocity  but  includes  a  capability  to  consider 
aerodynamic  lift  and  drag  as  well  as  a  nose  and  main  landing  gear  set 
with  weight  distribution  controlled  by  eg  location.  Finally  TAKOFF 
is  a  program  which  determines  the  Lakeoff  distance  required  on  soil 
and  on  pavement  for  an  aircraft  with  all  the  variables  of  the  first 
two  programs  but  requires  a  precise  description  of  thrust  as  a  func  in 
of  velocity.  The  program  uses  a  4-step  Runge-Kutta  integration  rout^..  t 
for  the  distance  computation.  All  the  programs  could  have  been  in¬ 
corporated  into  a  single  program,  however,  to  avoid  a  number  of  confus¬ 
ing  options,  separate  codes  were  used  for  each. 

Complete  documentation  for  use  of  the  three  programs  with  a  sample 
problem  for  each  is  given  in  the  appendices.  Appendix  A  contains  a 
complete  listing  of  the  three  codes.  Appendix  B  gives  instructions  for 
their  use  with  sample  problems  to  illustrate  the  set-up  of  the  input, 
data  decks.  Results  of  the  sample  problems  are  contained  in  Appendix  C 
with  a  plotted  comparison  of  results. 

To  assure  that  the  programs  were  working  properly,  the  same  simula¬ 
tion  which  was  made  by  Boeing  is  reproduced  as  closely  as  possible. 

This  correlation  as  well  as  estimates  of  the  performance  capability  of 
the  C-141  and  C-5  aircraft  are  presented  in  the  next  sections. 
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SECTION  IV 


CORRELATION  OF  COMPUTER  PROGRAM  RESULTS 
USING  THE  BOEING  367-80  AIRCRAFT 


1.  Discussion  -  To  assure  that  the  computer  programs  were  accurately 
coded,  a  performance  correlation  was  made  with  results  obtained  by 
Boeing  using  the  367-80  aircraft  as  the  analytical  vehicle.  Comparisons 
were  made  of  paved  takeoff  performance,  soil  takeoff  performance  and 
velocity  histories  of  soil  drag  ratio  and  wheel  sinkage. 

2.  Results  -  Figures  12,  13  and  14  show  paved  surface  performance  of 
the  367-80  operating  on  Boeing  Field.  The  analytical  data  agrees  well 
with  experimental  results  and  would  probably  come  even  closer  had  the 
Boeing  Field  runway  slope  been  accounted  for.  Figure  15  shows  analy¬ 
tical  and  experimental  drag  ratio  comparisons.  Drag  ratio  is  simply 
the  total  soil  drag  divided  by  the  total  vertical  load  (weight  minus 
lift).  The  difference  between  experimental  and  analytical  results 

in  this  figure  is  quite  sizable  at  some  speeds.  Discussions  with  Boeing, 
however,  indicate  that  the  experimental  data  were  quite  scattered  for 
this  test  and  that  under  the  circumstances,  the  analytical  results  are 
reasonable  approximations.  Figures  16  and  17  are,  respectively, 
analytical  sinkage  history  and  a  velocity-distance  comparison  for  the 
runs  of  Figure  15.  Figure  18  shows  a  drag  ratio  history  similar  to 
that  of  Figure  15  except  for  a  sandy  soil.  Figures  19  and  20  show 
the  percentage  increase  in  takeoff  distance  as  a  function  of  soil 
strength.  Figure  19  shows  a  considerable  variation  between  the  Boeing 
results  and  the  FDDS  correlation.  There  are  a  number  of  reasons  for 
this  variation,  but  the  greatest  cause  is  probably  differences  in  the 
data  used  for  the  analytical  vehicle.  Slight  variations  in  tire 
deflection  curves,  center  of  gravity  location  and  other  similar 
parameters  could  very  easily  cause  the  magnitude  of  the  variations 
shown.  The  data  for  the  367-80  vehicle  used  in  this  correlation  are 
presented  in  Appendix  D.  The  367-80  correlation  study  has  given 
sufficient  confidence  in  the  method  and  the  computer  programs  to 
warrant  making  a  study  of  the  C-141  and  (.-5  vehicles  which  will  be 
required  to  operate  on  substandard  landing  sites.  Assessment  of  the 
performance  capability  of  the  C-141  and  C-5  aircraft  are  covered 
in  detail  in  the  next  two  sections. 
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GROUND  ROLL  TAKEOFF  DISTANCE  (FT) 


120 


GROUND  DISTANCE  (FT) 


Figure  13  Paved  Surface  Ground  Velocity  as  a  Function  of  Ground 
Distance,  367-80,  180,300  lb.  Gross  Weight 
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GROUND  VELOCITY  (KNOTS) 


GROUND  DISTANCE  (FT) 


Figure  1A  Paved  Surface  Ground  Velocity  as  a  Function  of  Ground 
Distance,  367-80,  157,200  lb.  Gross  Weight 
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Figure  17  Cround  Velocity  as  a  Function  of  Ground  Distance  for 
367-80  at  138,700  lb.  Cross  Weight. 
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SECTION  V 


PERFORMANCE  OF  C-141  AIRCRAFT  OPERATING  ON 
CLAY  AND  SAND  AIRFIELDS 


1.  Discussion  -  Since  the  C-141  aircraft  is  one  of  the  most  widely 
used  cargo  and  troop  transport  vehicles  presently  being  used  by  the  Air 
Force,  it  follows  that  questions  continue  to  arise  concerning  its 
ability  to  operate  from  landing  fields  which  are  not  paved.  The  results 
presented  in  the  next  paragraph  attempt  to  make  an  assessment  of  the 
drag  loads  experienced  and  the  change  in  takeoff  performance  which  will 
result  from  operation  on  clay  and  sand  type  fields.  No  attempt  was 
made  to  determine  the  number  of  times  a  particular  takeoff  run  could 

be  made  on  a  single  soil  strength  before  the  surface  strength  character¬ 
istics  deteriorated  to  an  unacceptable  level.  This  study  is  concerned 
only  with  obtaining  a  minimum  soil  strength  for  a  particular  configura¬ 
tion.  The  basic  data  used  in  the  computer  program  is  given  in 
Appendix  E. 

2.  Results  -  As  in  the  previous  section,  the  results  of  operating  the 
C-141  aircraft  on  clay  and  sand  soil  is  compared  to  paved  surface  opera¬ 
tion,  and  therefore,  paved  surface  takeoff  performance  is  presented 
first  (Figure  21).  Exact  duplication  of  the  takeoff  distance  presented 
in  the  C-141  handbook  (Reference  5)  could  not  be  made  since  a  more 
sophisticated  outer  program  war  used  for  the  handbook  data.  The 
handbook  resulr  ire  calculated  assuming  that  a  rotation  velocity 

and  a  takeoff  velocity  must  be  reached  during  the  ground  roll.  This 
standard  technique  requires  a  time  history  of  angle  of  attack  and 
results  in  a  varying  lift  and  drag  coefficient  between  the  rotation 
speed  and  liftoff  speed.  The  present  program  assumed  that  no  angle 
of  attack  change  occurs  until  liftoff  speed  is  reached.  As  a  result, 
most  of  the  program  results  show  a  shorter  paved  surface  ground  distance. 

To  obtain  a  representative  sample  of  C-141  soil  performance,  calcula¬ 
tions  were  made  for  three  gross  weights,  three  soil  strengths  for  clay, 
three  soil  strengths  for  sand  and  tve  tire  pressures.  Table  1  lists 
the  figure  number  for  each  configuration  analyzed.  he  figu.es  are 
arranged  in  groups  o*'  three,  ’the  first  of  each  group  is  the  ratio  of 
soil  drag  and  vertical  load  as  a  Junction  of  ground  speed,  the  second  is 
main  gear  sinkage  for  the  same  condition  and  the  third  of  the  series  is 
ground  velocity  as  a  function  of  ground  distance.  At  the  end  of  each 
three  gross  weight  series  is  a  curve  of  percent  increa.se  in  takeoff 
distance  over  paved  surface  as  a  function  of  soil  strengths.  Figures 
34,  44,  *>4  and  64  are  the  four  such  curves  produced  in  .ills  study. 

Figure  34  and  Figure  44  may  be  compared  to  show  the  effects  cf  a  change 
in  tire  pressure  on  c!«-  and  Figures  54  and  64  for  the  same  results  on 
sand.  Comparison  of  rigures  34  and  54  will  give  the  results  of  changing 
from  clay  fo  sand  so:  1  for  the  sa~c  condition.  Although  ?  :e  present  study 
was  limited,  it  appears  that  j  reasonably  good  off-runway  capability 
exists,  particularly  in  the  higher  cone  index  strengths  <  r'vest  igated . 
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TABLE  1 


C-lUl  CONFIG  /RATION  SUMMARY 


Gross  Weight 

Surface  Type 

Surface  Strength 

Tire  Pressure 

Figure  Nr 

160,000 

Paved 

(Cl  for  Clay) 

(  G  for  Sand) 

2? 

220, OuO 

Paved 

- 

- 

23 

280,000 

Paved 

- 

- 

2k 

160,000 

Clay 

150,200,300 

100  Psi 

25 

160,000 

Clay 

150,200,300 

100  Psi 

26 

160,000 

Clay 

150,200,300 

100  Psi 

27 

220,000 

Clay 

220,260,325 

100  Psi 

28 

220,000 

Clay 

220,260,325 

100  Psi 

29 

220,000 

Clay 

220,260,325 

100  Psi 

30 

280,000 

Clay 

275,300,350 

100  Psi 

31 

280,000 

Clay 

275,300,350 

100  Psi 

32 

280,000 

Clay 

275,300,350 

100  Psi 

33 

160,000 

Clay 

150,200,300 

130  Psi 

35 

160,000 

Clay 

150,200,300 

130  Psi 

36 

160,000 

Clay 

150,200,300 

130  Psi 

37 

220,000 

Clay 

220,260,325 

130  Psi 

38 

220,000 

Clay 

220,260,325 

130  Psi 

39 

220,000 

Clay 

220,260,3  .’5 

130  Psi 

1*0 

260,000 

Clay- 

275, 300,3 50 

130  Psi 

1*1 

280,000 

Clay 

275,300,350 

130  Psi 

1*2 

280,000 

Clay 

275,300  350 

130  Psi 

U3 

160,000 

Sand 

U5.65.I00 

iOO  Psi 

U5 

160,000 

Sand 

1*5,65,100 

100  Psi 

1*6 

100,000 

Sand 

1*5,65,100 

100  Psi 

1*7 

220,000 

Sand 

60,80,100 

100  Psi 

1*8 

220,000 

OUT,  j. 

00, 80, 100 

100  Psi 

1*9 

220,000 

Sand 

60,80,100 

100  Psi 

50 

280,000 

Sand 

70,100,125 

100  Psi 

51 

280,000 

Sand 

70,100,125 

100  Psi 

52 

280,000 

Sand 

70,100,125 

IOO  Psi 

53 

IbO ,0 JO 

Sand 

1*5,65,100 

130  Psi 

55 

ioO ,000 

Sand 

1*5,65,100 

130  Psi 

56 

loo ,000 

Sand 

1*5,65,100 

130  Psi 

57 

2.  0,000 

Sand 

60 ,80 ,100 

130  Psi 

58 

220  ,000 

Sand 

60,80,100 

130  Psi 

59 

220 ,000 

Sand 

oO, 80, 100 

130  Psi 

60 

2 'iO  ,OUO 

Sand 

70,100,125 

130  Psi 

61 

280  ,000 

Sand 

70,100,125 

130  Psi 

62 

200,000 

Sand 

70,100,125 

130  Psi 

63 

35 


Figure  22.  Airplane  Velocity  as  a  Function  of  Ground 
Distance,  C-Hl,  160,000  lb.  Cross  Weight 
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Figure 


Airplane  Velocity  as  a  Function  of  Ground 
Distance,  C-141,  220,000  lb.  Crus j  Weiguc 
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Figure  24.  Airplane  Velocity  as  a  Function  of  Cround 
Distance,  C-141,  280,000  lb.  Cross  Weight 
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Figure  26  Soil  Sinkage  as  a  Function  of  G.ound  Velocity 
for  C-141  at  160,000  Lb.  Gross  Weight 
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Figure  27.  Airplnne  Velocity  as  a  Function  of  Cround 

Distance  for  C-141  at  160,000  lb.  Gross  Weight 
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Figure  28  Drag  Ratio  as  a  Function  of  Ground  Velocity  for 
C-’  tl  at  220,000  Lb.  Gross  Weight 
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Sinkage  as  a  Function  of  Ground  Velocity 
C-141  at  220,000  Lb.  Gross  Weight 
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Distance  for  C-141  at  780,000  lb.  Gross  Weight 


Figure  34  Percent  Increase  in  Takeoff  Distance  Over 
Paved  Surface  for  C-141  Aircraft 
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Figure  37.  Airplane  Velocity  as  a  Function  of  Ground  Distance 
for  C-141  at  160,000  lb.  Gross  Weight 
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Figure  38  Drag  Ratio  as  a  Function  of  Groun^  Velocity  for 
C-141  at  220,000  Lb.  Gross  Weight 
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Figure  40.  Airplane  Velocity  as  a  Function  of  Ground  Distance 
for  C-141  at  220,000  lb.  Gross  Weight 
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atio  as  a  Function  of  Ground  Velocity 
141  at  280,000  Lb.  Gross  Weight 
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jure  42  Soil  Sinkage  as  a  Function  of  Ground  Velocity 
for  C-I41  at  280,000  Lb.  Gross  Weight 
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Figure  44  Percent  Increase  in  Takeoff  Distance  Over  Paved 
Surface  for  C-141  Aircraft 
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Figure  47.  Airplane  Velocity  as  a  Function  of  Ground  Distance 
for  C-141  Aircraft  at  160,000  lb.  Gross  Weight 
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figure  50.  Airplane  Velocity  as  a  Function  of  Ground  Distance 
for  C-141  at  220,000  lb.  Gross  Weight 
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Figure  52  Soil  Sinkagc  as  a  Function  of  Ground  Velocity  for 
C-141  at  260,000  Lb.  Gross  Weight 


Figure  53.  Airplane  Velocity  as  a  Function  of  Ground  Distance 
for  C-141  at  280,000  lb.  Gross  Weight 
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Figure  54  Percent  Increase  in  Takeoff  Distance  Over  Paved 
Surface  for  C-141  Aircraft 


Figure  56  Soil  Sinkage  as  a  Function  of  Cround  Velocity 
for  C-141  at  160,000  Lb.  Gross  Weight 


Figure  57.  Airplane  Velocity  as  a  Function  of  Ground  Distance 
for  C-141  at  160,000  lb.  Gross  Weight 
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Figure  59  Soil  Sinkage  as  a  Function  of  Ground  Velocity 
for  C-141  at  220,000  Lb.  Gross  Weight 
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Figure  60.  Airplane  Velocity  as  a  Function  of  Ground  Distance 
for  C-141  at  220,000  lb.  Gross  Weight 
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Figure  62  Soil  Sinkage  as  a  Function  of  Ground  Velocity 
for  C-141  at  280,000  Lb.  Gross  Weight 
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Figure  63.  Airplane 
for  C-141 
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Figure  64  Percent  Increase  in  Takeoff  Distance  Over  Paved 
Surface  for  C-141  Aircraft 


SECTION  VI 


PERFORMANCE  OF  C-5  AIRCRAFT  OPERATING  ON  CLAY  AND  SANDY  AIRFIELDS 

Discussion  -  One  of  the  basic  design  requirements  of  the 
C-5  aircraft  is  to  operate  on  semi -prepared  soil  airfields.  To 
satisfy  this  requirement,  some  means  of  predicting  the  takeoff  distance 
analytically  is  necessary  for  construction  of  flight  performance  hand¬ 
book  charts,  since  full  scale  testing  for  each  configuration  is 
prohibited  from  a  cost  standpoint.  The  curves  presented  in  this 
section  are  of  the  type  which  could  provide  handbook  information. 

The  curves  are,  of  course,  based  only  on  a  limited  test  program  for 
full  scale  aircraft  but  reasonably  accurate  predictions  appear  to  be 
possible. 


Results  -  The  format  for  presenting  the  C-5  data  is  the  same 
as  that  for  the  C-141.  A  paved  surface  takeoff  distance  curve  is 
presented  in  Figure  65  for  a  reference  since  all  soil  takeoff  distance 
curves  (Figures  78,  88,  98  and  108)  are  presented  as  a  percent  Increase 
in  the  paved  distance.  Paved  takeoff  distances  calculated  with  the 
program  were  shorter  than  the  values  calculated  by  Lockheed  because 
of  our  simplification  of  the  rotation  phase  by  elimination  of  the 
varying  lift  and  drag  coefficients  (see  Appendix  F  for  further  details 
and  computer  data) . 


Table  II  gives  an  index  of  the  data  presented 
for  the  C-5.  Two  soil  types  with  two  tire  pressures,  three  gross 
weights  and  three  soil  strengths  were  used  to  obtain  the  figures. 
Summary  figures  (78,  88,  98  and  108)  may  be  compered  to  see  the 
effects  of  parameter  change  on  the  increase  in  takeoff  distance. 
Figures  78  and  88  show  the  difference  in  takeoff  distance  because 
of  a  change  in  tire  pressure  on  clay.  Figures  98  and  108  show  a 
similar  change  for  sand.  Comparison  of  Figures  78  and  98  will  show 
the  effect  of  changing  from  clay  to  sand  for  the  same  weight  and  tire 
pressure  conditions.  Figures  88  and  108  give  a  comparison  of  changing 
from  clay  to  sand  tor  a  higher  tire  pressure. 
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TABLE  2 


C-5  CONFIGURATION  SUMMARY 


Gross 

Weight 

Surface 

Type . . . 

Surface 

Strength 

Tire 

Pressure 

Figure 

Number 

400,000 

Lb. 

Paved 

(Cl  for  Clay) 

(  G  for  Sand) 

66 

500,000 

Lb. 

Paved 

- 

- 

67 

571,000 

Lb. 

Paved 

- 

- 

68 

400,000 

Lb. 

Clay 

175,200,250 

95  Psi 

69 

400,000 

Lb. 

Clay 

175,200,250 

95  Psi 

70 

400,000 

Lb. 

Clay 

175,200,250 

95  Psi 

71 

500,000 

Lb. 

Clay 

200,225,275 

95  Psi 

72 

500,000 

Lb. 

Clay 

200,225,275 

95  Psi 

73 

500,000 

Lb. 

Clay 

200,225,275 

95  Psi 

74 

571,000 

Lb. 

Clay 

225,250,300 

95  Pai 

75 

571,000 

Lb. 

Clay 

225,250,300 

95  Psi 

76 

571,000 

Lb. 

Clay 

225,250,300 

95  Psi 

77 

400,000 

Lb. 

Clay 

175,200,250 

120  Psi 

79 

400,000 

Lb. 

Clay 

175,200,250 

120  Psi 

80 

400,000 

Lb. 

Clay 

175,200,250 

120  Psi 

81 

500,000 

Lb. 

Clay 

200,225,275 

120  Pai 

82 

500,000 

Lb. 

Clay 

200,225,275 

120  Psi 

83 

500,000 

Lb. 

Clay 

200,225,275 

120  Psi 

84 

571,000 

Lb. 

Clay 

225,250,300 

120  Psi 

85 

571,000 

Lb. 

Clay 

225,250,300 

120  Psi 

86 

571,000 

Lb. 

Clay 

225,250,300 

120  Psi 

87 

400,000 

Lb. 

Sand 

75,100,125 

95  Psi 

89 

400,000 

Lb. 

Sand 

75,100,125 

95  Psi 

90 

400,000 

Lb. 

Sand 

75,100,125 

95  Psi 

91 

500,000 

Lb. 

Sand 

80,110,150 

95  Psi 

92 

500,000 

Lb. 

Sand 

80,110,150 

95  Psi 

93 

500,000 

Lb. 

Sand 

80,110,150 

95  Psi 

94 

571,000 

Lb. 

Sand 

85,120,175 

95  Psi 

95 

571,000 

Lb. 

Sand 

85,120,175 

95  Psi 

96 

571,000 

Lb. 

Sand 

85,120,175 

95  Psi 

97 

400,000 

Lb. 

Sand 

75,100,125 

120  Pai 

99 

400,000 

Lb. 

Sand 

75,100,125 

120  Psi 

100 

400,000 

Lb. 

Sand 

75,100,125 

120  Psi 

101 

500.000 

Lb. 

Sand 

80,110,150 

120  Pel 

102 

500,000 

Lb. 

Sand 

80,110,150 

120  Psi 

103 

500,000 

Lb. 

Sand 

80,110,150 

120  Psi 

104 

571,000 

Lb. 

Sand 

85,120,175 

120  Psi 

105 

571,000 

Lb. 

Sand 

85,120,175 

120  Pr.i 

106 

571,000 

Lb. 

Sand 

85,120,175 

120  Psi 

107 

80 
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Figure  66.  Airplane  Velocity  as  i  Function  of  Ground  Distance 
for  C-5  at  400,000  lb.  Gross  Weight 
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Figure  67.  Airplane  Velocity  as  a  Function  of  Ground  Distance 
for  C-5  at  500,000  lb.  Gross  Weight 
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Figure  68.  Airplane  Velocity  as  a  Function  of  Ground  Distance 
for  C-5  at  571,000  lb.  Gross  Weight 
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Figure  72  Drag  Ratio  as  a  Function  of  Ground  Velocity  for 
C-5  at  500,000  Lb.  Gross  Weight 


Airplane  Velocity  as  a  Function  of  Ground  Distance 
for  C-5  at  500,000  lb.  Gross  Weight 


o 

o 


w 

o 


H 

(A 


O 

« 

o 


9  j 


> 


Ftaure  77  •  Airplane  Velocity  as  a  Function  of  Ground  Distance 
for  c-5  at  571,000  lb.  Gross  Weight 
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Figure  78  Percent  Increase  in  Takeoff  Distance  Over  Paved 
Surface  for  C-5  Aircraft 
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Figure  83  Soil  Sinkage  as  a  Function  of  Ground  Velocity 
for  C-5  at  500,000  Lb.  Gross  Weight 


100 


000  Lb.  Gros-s  Weight 
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Figure  86  Soil  Sinkage  as  a  Function  of  Ground  Velocity 
for  C-5  at  571,000  Lb.  Gross  Weight 
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at  571,000  lb.  Gross  Welgl 
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Figure  88  Percent  Increase  in  Takeoff  Distance  Over  Paved 
Surface  for  C-5  Aircraft 


Figure  89  Drag  Ratio  as  a  Function  ->f  Ground 
for  C-5  at  400,000  Lb.  Gross  Weight 
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GROUND  DISTANCE  (100  FT) 


Figure  91.  Airplane  Velocity  as  a  Function  of  Ground  Distance 
for  C-5  at  400,000  lb.  Gross  Weight 
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Figure  92  brag  Ratio  as  a  Function  of  Ground  Velocity  for 
C-.'5  at  500,000  Lb.  Gross  Weight 
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Figure  98  Percent  Increase  in  Takeoff  Distance  Over  Paved 
Surface  for  C-5  Aircraft 


(dVOl  'lV)UH3A/DYHa  HY3D)  OIXVH  OVttd 


115 


>i 


Figure  99  Drag  Ratio  as  a  Function  of  Ground  Velocity  for 
C-5  at  400,000  Lb.  Gross  Weight 
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Figure  100  Soil  Sinkage  as  a  Function  of  Ground  Velocity  for 
C-5  at  400,000  Lb.  Gross  Weight 
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C-5  at  400,000  lb 
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Figure  102  Drag  Ratio  as  a  Function  of  Ground  Velocity  for 
C-5  at  500,000  Lb.  Gross  Weight 


Figure  103  Soil  Sinkage  as  a  Function  of  Ground  Velocity 
for  C-5  at  500,000  Lb.  Gross  Weight 


Figure  104.  Airplane  Velocity  as  a  Function  of  Ground  Distance 
for  C-*  at  50U.UU0  lb.  Gross  Weight 
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SECTION  VII 


AREAS  FOR  FURTHER  DEVELOPMENT 

Discussion  -  In  a  technology  such  as  the  one  described  in  the 
preceding  sections,  there  always  exists  a  number  of  areas  where 
because  of  limited  time  or  lack  of  better  scientific  evidence, 
assumptions  and  simplifications  are  made  in  the  techniques  used. 

In  the  case  of  predicting  soil  drag  and  soil  sinkage  under  dynamic 
deformation,  a  method  was  developed  based  on  empirical  curves 
which  forced  analytical  results  to  conform  to  a  limited  amount  of 
test  data.  In  this  section  we  will  attempt  to  point  out  a  few  areas 
whore  improvement  and  extended  verification  of  techniques  and 
assumptions  should  be  made. 

Areas  Where  Further  Development  is  Required  -  A  number  of  areas 
have  been  observed  where  further  study  could  provide  a  more  complete 
understanding  of  the  mechanics  of  dynamic  soil  deformation  caused  by 
a  rolling  tire.  These  areas  with  some  suggested  avenues  of  improvement 
are  discussed  below. 

2.1  Mobility  Number  -  The  present  method  uses  the  concept  of 
Mobility  Number  developed  by  WES  to  predict  soil  sinkage  and  drag. 

The  results  of  plotting  mobility  number  for  various  tire,  load  and 

soil  strength  combinations  has  shown  a  great  deal  of  scatter.  The  tires 
used  were  for  both  aircraft  and  land  vehicles.  An  effort  to  minimize 
this  mobility  number  scatter  by  considering  only  aircraft  tires  for 
example,  could  significantly  improve  the  sinkage  and  drag  predictions. 

2.2  Dynamic  Factor  -  The  dynamic  factor  is  a  parameter  which 
has  been  established  to  correct  for  high  speed  effects  when  scaling 

drag  and  sinkage  values  based  on  low  speed  tests  to  the  range  of  interest 
for  aircraft.  The  curve  of  Figure  6  was  drawn  based  on  a  limited 
amount  of  test  data  most  of  which  was  in  the  low  speed  range.  A  simple 
study  of  some  preliminary  data  obtained  (Reference  6)  has  shown 
considerable  deviation  from  Figure  6,  particularly  at  high  .speeds. 
Additional  studies  of  this  dynamic  deformation  effect  may  prove  fruitful 
in  improving  prediction  accuracy. 

Tire  Lift  Curve  -  Toe  assumption  has  been  made  that  the  soil 
exerts  a  dynamic  pressure  and  a  corresponding  lift  on  the  tire  much  iike 
an  aircraft  wing.  Consequently,  values  of  an  associated  tire  lift 
coefficient  were  determined  based  on  a  point  at  which  the  tire  reaches 
the  planing  velocity,  and  significant  reduction  in  drag  is  observed. 

The  value  of  lift  coefficient  controls  the  peak  of  the  drag  ratio 
velocity  histories  presented  in  previous  sections  and,  as  such,  becomes 
quite  important  in  the  overall  prediction  technique.  More  studies  of 
th's  "planing"  phenomena  and  the  variation  of  lift  coefficient  with  tire 
configuration  are  required  before  a  clear  unuet standing  of  the  effects  can 
be  attained. 
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2.4  Impingement  Drag  Coefficient  -  A  parameter  closely 
associated  with  the  lift  coefficient  is  impingement  drag  coefficient. 
This  parameter  shown  in  Figure  11  is  controlled  by  lift  coefficient 
through  the  planing  velocity.  The  Impingement  drag  coefficient  has 
a  significant  Influence  on  the  shape  of  the  drag  ratio-velocity  curves. 
The  assumption  of  a  constant  level  up  to  planing  velocity  and  a  low  but 
constant  value  beyond  1.5  times  the  planing  velocity  may  be  too  severe. 
Future  studies  will  refine  the  curve  shape  to  conform  to  actual  test 
results. 


2.5  Soil  Strength  Measurements  -  One  parameter  which  affects  the 
overall  technique  is  the  soil  strength.  A  method  is  needed  which  will 
provide  a  simple  repeatable  measurement  of  soil  strength  for  all 
types  of  uoll.  Present  methods  are  either  cumbersome  or  operator- 
dependent  and,  as  sucn,  provide  a  major  hindrance  to  future  soil 
dynamics  studies. 
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SECTION  VIII 


CONCLUSIONS 


1.  Computer  Programs  -  Three  computer  programs  have  been  developed 
which  will  predict  soil  slnkage  and  resulting  drag  loads  for  an  aircraft 
operating  on  clay  and  sandy  type  airfields.  The  three  programs  will 
allow  assessment  of  loads  and  performance  for  varying  quantities  of 
vehicle  descriptive  data.  The  data  requirements  specified  in  this 
fashion  will  permit  assessment  of  a  wide  variety  of  vehicles  and  con¬ 
figurations  and  permit  a  "quick  look"  or  "detailed  assessment"  depend¬ 
ing  on  the  amount  of  time  available  to  determine  the  necessary 
descriptive  data. 

The  program  appear  to  give  results  which  compare  favorably  with 
the  limited  exper iz.ental  data,  although  further  refinement  and  accuracy 
improvement  v'll  undoubtedly  result  when  findings  of  a  research  program 
presently  underway  are  incorporated . 

Additional  consideration  should  be  given  to  the  following  areas 
which  are  the  basic  weaknesses  in  the  generality  of  the  programs. 

a.  Zero  initial  velocity  cannot  be  used  in  the  TAKOFF  program 
because  it  results  in  a  divergent  solution  caused  by  high  zero  velocity 
drag  values  which  override  thrust.  Studies  of  initial  loads  to  begin 
motion  for  a  static  soil  slnkage  case  would  provide  an  insight  into 
this  problem.  Beginning  calculation  with  two  or  three  ft/sec  initial 
velocity  sr  s  to  solve  this  problem  with  TAKOFF. 

b.  The  -grams  assume  that  regardless  of  the  soil  strength  the 
minimum  soil  .  'cage  is  0.1  inch.  This  lower  limit  has  been  established 
to  avoid  negatx  slnkage  coefficients  as  vertical  loads  becomes 
smaller  near  the  end  of  a  takeoff  roll.  As  can  be  seen  from  Figures 

34,  78  and  other  increased  takeoff  distance  curves,  the  percent 
increase  approaches  3  to  5  percent  at  high  soil  strengths.  The  fixed 
slnkage  limit  is  partially  responsible  for  the  increased  takeoff 
distance  remaining  at  3  to  3  percent  even  at  very  high  soil  strengths 
where  normally  one  would  expect  it  to  approach  zero.  The  limit 
could  be  lowered  slightly  but  with  questionable  benefit  since  other 
factors  such  as  a  minimum  rolling  resistance  coefficient  of  0.03  also 
has  an  effect.  This  level  was  established  on  the  basis  that  soil 
rolling  resistance  regardless  of  soil  strength  could  not  be  less  than 
the  rigid  surface  rolling  resistance.  Additional  testing  may  clarify 
this  situation  somewhat. 

c.  Another  area  whore  difficulty  is  encountered  in  using  the  TAKOFF 
program  is  in  determining  the  range  of  soil  strengths  where  realistic 
results  can  be  obtained.  If,  for  example,  too  low  a  value  of  soil 
strength  is  used  the  solution  will  diverge  or  reach  some  steady  state 
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condition.  Some  consideration  was  given  to  the  solution  to  this 
problem  but  no  changes  were  aade  in  the  program  to  correct  it.  The 
method  used  in  the  present  analysis  was  to  begin  the  calculation  with 
a  very  high  soil  strength  such  as  400  Cl  or  200  G  and  work  to  the  lower 
strengths  by  checking  the  percent  increase  in  takeoff  diatance  over 
paved  surface.  Percent  Increase  values  approaching  100  we- e  considered 
to  give  a  practical  lower  limit  on  soil  strength  for  the  program.  This 
has  worked  satisfactorily  sine?,  no  problem  occurs  until  a  lower  limit 
of  approximately  80  Cl  or  20  G  is  tried. 

2.  Method  -  The  method,  although  in  need  of  continued  refinement,  appears 
to  give  reasonably  accurate  results  end,  as  such,  provides  the  user  with 
a  technique  for  evaluating  a  number  of  technical  problems  associated 
with  ground  operation  of  aircraft  which  were  heretofore  only  estimated 
with  engineering  judgment.  Some  areas  where  fruitful  results  can  be 
obtained  with  the  method  are  listed  below. 

a.  Determine  takeoff  performance  on  soil  airfields  for  a  particular 
configuration. 

b.  Determine  design  load  data  for  landing  gear  expected  to  be  used 
on  soil  airfields. 

c.  Optimize  tire  sizes  for  a  particular  aircraft  design. 

d.  Optimize  tire  geometry  and  number. 

e.  Calculate  dynamic  effects  of  operation  on  soil  airfields. 

It  is  expected  therefore  that  in  the  future,  techniques 
such  as  the  one  given  herein  will  be  used  as  a  standard  design 
calculation  for  aircraft  required  to  operate  on  improvised  airfields. 
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COMPUTER  PROGRAM 
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APPENDIX  B 


SAMPLE  INPUT  FOR  L.-IL  DRAG  PROGRAMS 


1.  Discussion  -  This  Appendix  gives  the  input  dais  setup  for  a  samp  e 
problem  wh'-r  'll  be  solved  using  each  of  the  Lhree  programs.  The 
problem  is  to  determine  the  performance  capability  of  the  Boeing  367  80 
aircraft  at  150,000  lb.  gross  weight  operating  on  clay  soil  with  a  Cone 
Index  of  90.  The  problem  will  be  solved  assuming  three  levels  of  data 
compl.- 1- ness .  The  first  (Problem  '.)  assumes  that  all  data  given  in 
Appendix  C  is  known  and  that  the  takeoff  distance  time  history  on  clay 
and  on  paved  surface  is  required  in  order  to  determine  the  increase  in 
takeoff  distance  which  will  result  from  operation  on  the  clay  soil. 

The  second  problem  (Probler.  B)  will  assume  that  no  engine  thrust  data 
is  available.  This,  of  course,  precludes  the  time  history  and  the 
-associated  takeoff  distance  but  will  give  the  sane  soil  sinkage  ar.  1 
drag  which  was  obtained  from  the  more  complete  problem.  The  final 
case  (Problem  C)  assumes  that  only  the  gross  weight  of  the  vehicle  and 
the  number  ard  type  of  tires  used  on  the  vehicle  are  available.  This 
solution  neglects  aerodynamic  effects,  and  therefore  results  in  slightly 
higher  drag  loads  tl  .-n  the  other  program  but  is  conservative  for 
design  and  has  the  advantage  of  requiring  a  minimum  of  input  data. 

Data  and  deck  setup  for  each  problem  will  now  be  covered  as  a  separate 
paragraph. 

2.  Problem  A  -  Using  the  data  cf  Appendix  D  for  the  367-80  aircraft 
and  the  computer  program  TAKOFF  given  in  Appendix  A,  calculate  a 
tskeoff  time  history  for  tu  ‘  aircraft  at  150,000  lb.  gross  weight 
operating  on  a  clay  soil  with  a  rone  index  of  90,  and  a  takeoff  time 
history  on  a  paved  runway. 

The  input  data  consists  of  seven  distinct  sections  as  shown  in 
Figure  109  The  details  cf  each  section  are  as  follows: 


Section  1  (one  card)  -  A  rodynamic  Pat? 


Card 

Column 

Format 

Variable 

Name 

Def inition 

1-12 

F12.4 

AW 

Wing  Area  (Ft2) 

13-24 

F12.4 

CL 

Takeoff  Lift  Coefficient 

25-36 

F12.4 

CD 

Takeoff  Drag  Coefficient 

37-48 

F12.4 

SIG 

Atmospheric  Density  Patio 

49-60 

112 

K 

Number  of  Values  Each  of  Thrusr 

and  Velocity 
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Section  2  -  Thrust  Data 


Card 

Variable 

Column 

Format 

Name 

Definition 

1-72 

6F12.4 

XVEL (I) 

K  Values  of  Velocity  (Knots) 

1-72 

6F12.4 

YTH(l) 

K  Values  of  Single  Engine  Thrust 

Corresponding  to  XVEL  (1)  Values 


Figure  109  Input  Data  Deck  for  Program  TAKOFF 
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Card  1 


Section  3  (2  Cards)  -  Wheel  and  Runway  Data 


Card 

Column 

Format 

Variable 

Name 

Definition 

1-12 

F12.4 

DWB 

Distance  Between  Nose  and  Main 

13-24 

F12.4 

WNG 

Gear  Struts  (Ft) 

Number  of  Nose  Wheels  Per  Strut 

25-36 

F12.4 

WMG 

(one  strut  allowed) 

Number  of  Main  Gear  Wheels  per 

37-48 

F12.4 

BN 

Side 

Width  of  Nose  Gear  Tire  (In) 

49-60 

F12.4 

DN 

Diameter  of  Nose  Gear  Tire  (In) 

61-72 

F12.4 

BM 

Width  of  Main  Gear  Tire  (In) 

Card  2 

Card 

Column 

Format 

Variable 

Name 

Definition 

1-12 

F12.4 

DM 

Diameter  of  Main  Gear  (In) 

13-24 

F12.4 

HTN 

Nose  Gear  Tire  Section  Height  (In) 

25-36 

F12.4 

HTM 

Main  Gear  Tire  Section  Height  (In) 

37-48 

F12.4 

RHOS 

Soil  (usually  2.6)  Density 

49-60 

F12.4 

SLOPE 

(Lb/Sec2/Ft4) 

Runway  Slope  (Rad) 

Section 

4  (1  Card) 

-  Integer  Control  Data 

Card 

Column 

Format 

Variable 

Name 

Definition 

1-12 

112 

NTANDN 

Number  of  Nose  Gear  Tandem  Wheel 

13-24 

112 

NTANDM 

Sets 

Number  of  Main  Gear  Tandem  Wheel 

25-36 

112 

NTYP 

Sets  Per  Side 

Soil  Type  (-1  for  Clay,  -2  for  Sand) 

37-48 

112 

NENGS 

Number  of  Engines 

49-60 

112 

N 

Number  of  Pairs  of  Load  and 

Deflection  Values  to  be  Read  for 
Nose  and  Main  Gear  Tire  Deflection 
Curves 
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Section  5  -  Tire  Deflection  Curves 


Card 

Column 

Format 

Variable 

Name 

Definition 

1-72 

6F12.4 

XFVN(I) 

N 

Values  of  Nose  Gear  Vertical 
Load  (Lb) 

1-72 

6F12.4 

YDFN(I) 

N 

Values  of  Nose  Gear  Tire 
Deflection  Corresponding  to 
the  XFVN(I)  Load  (In) 

1-72 

6F12.4 

XFVM(I) 

N  Values  of  Main  Gear  Vertical 
Load  (Lb) 

1-72 

6F12.4 

YDFM(I) 

N 

Values  of  Main  Gear  Tire 

Deflection  Corresponding  to 
the  XFVM(l)  Loads  (In) 


Section  6  -  Initial  Conditions  Data 


Card 

Column 

Format 

Variable 

Name 

Definition 

1-12 

F12.4 

SO 

Initial  Distance  for  Takeoff 
Integration  (Ft)  Usually  =  0 

13-24 

F12.4 

VO 

Initial  Velocity  for  Takeoff 
Integration  (Ft/Sec)  (Usually 

0  for  Paved  Runs  and  2  or  3 
for  Soil  Runs) 

25-36 

F12.4 

AO 

Initial  Acceleration  for  Takeoff 
Integration  (Ft/Sec^) 

Usually  ■  0 

37-48 

F12.4 

DT 

Integration  Interval  (Sec) 
(Usually  0.1  or  0.01) 

49-60 

F12.4 

ZS 

Initial  Soil  Sinkage  Estimate 

(In)  (Normal  Value  -  0.5) 
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Section  7  -  Run  Data 


Card 

Column 

Format 

Variable 

Name 

Definition 

1-12 

F12.4 

Cl 

Cone  Index  for  Clay  Soil  (PSI) 

Can  be  Zero  if  Paved  or  Sand 

Runs  are  Being  Made 

13-24 

F12.4 

G 

Cone.  Index  for  Sand  (PSI /In)  Can 
be  Zero  if  Paved  br  Clay  Runs 
are  being  Made 

25-36 

F12.4 

RRC 

Paved  Runway  Friction  Coefficient. 
Must  be  Zero  if  Clay  or  Sand 
runs  are  desired 

37-48 

F12.4 

W 

Total  Gross  Weight  of  Vehicle  (Lb; 

49-60 

F12.4 

DCG 

CG  Location  (Ft  Forward  of 

Main  Gear  Strut) 

61-72 

F12.4 

VT0 

Velocity  Required  for  Takeoff 
(Ft/Sec) 

For  additional  time  history  runs  on  the  same  soil  type,  only  the 
run  data  card  need  be  changed  to  reflect  a  new  soil  strength  or  gross 
weight  condition.  If  the  soil  type  is  to  be  changed,  however,  the 
Section  4  card  must  be  changed  to  reflect  the  desired  NTYP.  Tables 
3  and  4  give  a  listing  of  the  data  used  for  the  paged  surface  runs 
and  the  clay  soil  runs  respectively.  Notice  that  the  only  difference 
is  that  the  RRC  value  has  been  given  the  non-zero  value  of  0.02  and 
the  initial  velocity  has  been  set  to  zero  to  obtain  the  paved  takeoff 
time  history.  The  results  of  this  problem  are  given  in  Appendix  C. 

3.  Problem  B  -  Using  the  data  of  Appendix  B  for  the  367-80  aircraft 

and  the  computer  program  S0LDG2  given  in  Appendix  A,  compute  the  soil 

sinkage  and  drag  ratio  (total  drag  divided  by  total  vertical  load  -  lift  < 

minus  weight)  for  a  number  of  taxi  velocities  assuming  sufficient  thrust 

to  maintain  constant  speed.  The  aircraft  is  operated  on  clay  soil  with 

cone  index  of  90  and  gross  weight  of  150,000  lbs;. 

The  input  data  consists  of  four  sections  as  shown  in  Figure  110. 

Details  of  each  parameter  are  given  below: 


j 
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Section  1 

(3  Cards)  - 

Aircraft  and  Tire  Data 

Card  1 

Card 

Column 

Format 

Variable 

Name 

Definition 

1-12 

F12.4 

W 

Vehicle  Weight  (Lb) 

13-24 

F12.4 

ZS 

Initial  Estimate  of  Soil  Sinkage 
(In)  (Usually  0.5) 

25-36 

F12.4 

BN 

Diameter  of  Nose  Gear  Tire  (In) 

37-48 

F12.4 

HTN 

Nose  Tire  Section  Height  (In) 

49-60 

112 

N 

Number  of  Pairs  of  Nose  and  Main 
Tire  Deflection  Curves  to  be 
Read 

Figure  110  input  Data  Deck  for  Program  SOLDG2  and  SOLDG 


Card  2 


Card 

Variable 

Column 

Format 

Name 

Definition 

1-12 

F12.4 

AW 

Wing  Area  (Ft^) 

13-24 

F12.4 

BM 

Width  of  Main  Gear  Tire  (In) 

25-36 

F12.4 

DM 

Diameter  of  Main  Gear  Tire  (In) 

37-48 

F12.4 

HIM 

Main  Gear  Tire  Section  Height  (In) 

49-60 

F12.4 

CL 

Takeoff  Lift  Coefficient 

61-72 

F12.4 

CD 

Takeoff  Drag  Coefficient 

Card  3 

Card 

Variable 

Column 

Format 

Name 

Definition 

1-12 

F12.4 

DWB 

Distance  Between  Nose  and  Main 

Gear  Struts  (Ft) 

13-24 

F12.4 

DCG 

CG  Location  (Ft  Forward  of  Main 

Gear  Strut) 

25-36 

F14.4 

WNG 

Number  of  Nose  Wheels  Per  Strut 

(one  strut  allowed) 

37-48 

F12.4 

WMG 

Number  of  Main  Gear  Wheels  Per  Side 

49-60 

F12.4 

RHOS 

Soil  (usually  2.6)  Density 
(lb/secvft^) 

61-72 

F12.4 

SIG 

Atmospheric  Density  Ratio 

Section 

2  (1  Card) 

-  Control  Data 

Card 

Variable 

Column 

Format 

Name 

Definition 

1-12 

112 

NTYP 

Soil  Type  (-  1  for  Clay,  -2  for  Sand) 

13-24 

1.12 

NTANDN 

Number  of  Nose  Gear  Tandem  Wheel  Sets 

25-36 

112 

NTANDM 

Number  of  Main  Gear  Tandem  Wheel  Sets 

Per  Side 
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Section  3  Tire  Deflection  Data 


Card 

Column 

Format 

Var;’  e 

Na, 

Def inition 

1-72 

6F12.4 

XFVN(I) 

N 

Values  of  Nose  Tire  Vertical 
Load  (Lb) 

1-72 

6F12.4 

YDFN(I) 

N 

Values  of  Nose  Tire  Deflection 
Corresponding  to  Load  XFVN(I) 
(In) 

1-72 

6F12.4 

XFVM'I) 

N 

Values  of  Main  Tire  Vertical 
Load  (Lb) 

1-72 

6F12.4 

YDFM(I) 

N 

Values  of  Main  Tire  Deflection 
Corresponding  to  Load  XFVM(I) 
(In) 

Section  4  (1  Card 

Per 

Run)  Run  Data 

Card 

Column 

Format 

Variable 

Name 

Definition 

1-12 

F12.4 

CII 

Initial  Value  of  Soil  Strength 
Either  Cl  or  G  Depending  on 
Value  of  NTYP 

13-24 

F12.4 

DC  I 

Incremental  Soil  Strength 

25-36 

F12.4 

CIMAX 

Maximum  Value  of  Soil  Strength 
Desired 

37-43 

F12.4 

VI 

Initial  Taxi  Velocity  (Knots) 

49-60 

F12.4 

DELV 

Incremental  Taxi  Velocity  (Knots) 

61-72 

F12.4 

VMAX 

Maximum  Taxi  Velocity  Desired 

Qtnots) 


Using  the  run  data  of  Section  4,  the  program  will  determine  drag  ratio 
and  sinkage  for  each  taxi  velocity  for  a  particular  soil  strength  and 
will  then  increment  soil  strength  until  all  soil  strengths  have  been 
processed.  Table  5  gives  a  listing  of  the  data  used  in  the  problem. 
Results  &re  contained  in  Appendix  C. 


4.  Problem  C  -  Using  only  the  tire,  vehicle  weight  and  soil  data  from 
Appendix  D,  calculate  the  soil  sinkage  and  drag  ratio  values  corres¬ 
ponding  to  a  range  of  constant  taxi  velocities.  Use  computer  program 
SOLDC  from  Appendix  A  and  assume  clay  soil  with  cone  index  of  90  and 
vehicle  weight  of  150,000  lbs. 
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The  input  data  for  SOLDG  consists  of  four  sections  as  shown  in 
Figure  110.  Card  format  for  each  parameter  is  shown  below: 


Section  1 

(1  Card)  -  Vehicle  and  Tire  Data 

Card 

Column 

Format 

Variable 

Name 

Definition 

1-12 

F12.4 

W 

Wheel  Load  for  One  Wheel  of  One 
Main  Gear  Strut  (Lb) 

13-24 

F12.4 

ZS 

Initial  Estimate  of  Soil  Sinkage 
(In)  (Usually  0.5) 

25-36 

F12.4 

B 

Main  Gear  Tire  Width  (In) 

37-48 

F12.4 

D 

Main  Gear  Tire  Diameter  (In) 

49-60 

F12.4 

HT 

Main  Gear  Tire  Section  Height  (In) 

61-72 

F12.4 

RHOS 

Soil  Density  (Lb/Sec2/Ft4) 

Section  2  (1  Card) 

-  Control  Data 

Card 

Column 

Format 

Variable 

Name 

Definition 

1-12 

112 

N 

Number  of  Pairs  of  Main  Gear  Tire 
Deflection  Values  to  be  Read 

13-24 

112 

NTYP 

Soil  Type  (  «  1  if  Clay,  ■  2  if 
Sand) 

25-36 

112 

NTAND 

Number  of  Main  Gear  Tandem  Wheel 

Sets  Per  Strut 

Section  3  -  Tire 

Deflection  Data 

Card 

Column 

Format 

Variable 

Name 

Definition 

1-72 

6F12.4 

XFVN(I) 

N  Values  of  Main  Tire  Load  (Lb) 

1-72 

6F12.4 

YDFN(I) 

N  Values  of  Tire  Deflection  Which 

Correspond  to  the  XFVN(I)  Load 
(In) 


A 
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Card 

Column 

Format 

Variable 

Name 

Definition 

1-12 

F12.4 

CII 

Initial  Value  of  Soil  Strength 
Either  Clay  (Cl)  or  Sand  (G) 

13-24 

F12.4 

DC  I 

Incronental  Soil  Strength 

25-36 

F12.4 

CIMAX 

Maximum  Value  of  Soil  Strength 
Desired 

37-48 

F12.4 

VI 

Initial  Taxi  Velocity  (knots) 

49-60 

F12.4 

DELV 

Incremental  Taxi  Velocity  (Knots) 

61-72 

F12.4 

VMAX 

Maximum  Value  of  Velocity 

Desired  (Knots) 

Section  4  data  will  cause  the  program  to  calculate  soil  drag  and 
sinkage  for  various  velc  ;  let,  fnr  a  constant  ,joil  strength  and  then 
increment  soil  strength  ui-  il  all  strength  values  have  been  processed. 
Table  6  lists  the  input  data  for  Problem  C.  The  solution  to 
Problem  C  ie  given  in  Appendix  C. 
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TABLE  3 


SAMPLE  INPUT  DATA  FOR  PROBLEM  A  -  CLAY  AIRFIELD 


2621. 

0.641 

-20. 

0. 

40. 

70. 

9000. 

9000. 

13750. 

14000 

44. 

4. 

45. 

6.6 

0 

0. 

1250. 

0. 

.6 

0. 

2500. 

0.0 

1.2 

0. 

3.0 

90. 

0. 

0.048 

1. 

4. 

10. 

150. 

11000. 

12C00. 

14000. 

8. 

10.65 

10.75 

2.6 

2500. 

l 

3750. 

1.2 

1.7 

5000. 

10000. 

2.15 

3.75 

0. 

.1 

0. 

1 5000n 

9 

20. 

30. 

13000. 

13500 

30.7 

15.0 

0. 

5 

7500. 

3.7 

20000. 

6.9 

•  J 

2.36  219.57 


TABLE  4 

SAMPLE  INPUT  DAT A  FOR  PRLBLEM  A  -  PAVED  AIRFIELD 


2821. 

0.64  1 

-20. 

0. 

40. 

70. 

9C00 . 

9000. 

13750. 

14000 

44. 

4. 

45. 

6.6 

0 

0. 

1250. 

> 

s J  % 

.& 

0. 

2500. 

0.0 

1  .2 

0  . 

0.0 

O  • 

0. 

0.046 

i . 

4. 

10. 

150. 

UOCG. 

1 2  COO . 

ueco. 

8  . 

10.65 

10.75 

2.6 

2  400. 

1 

3  7  50. 

1.2 

1.7 

50  00. 

10000 

2.15 

3.  75 

0 . 

.  1 

o.\: 

i5c:o 

9 


2  C . 

30. 

13000. 

1  3500 

30.7 

0. 

15.0 

4 

7  >0?. 

3.  7 

2.000. 

6.9 
•  5 

5 

2.36 

219.5 
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APPENDIX  C 


SAMPLE  OUTPUT  FROM  SOIL  DRAG  PROGRAM 


1.  Discussion  -  Results  from  the  three  problems  presented  in  Appendix 
B  are  shown  as  Tables  7,  8,  9  and  10.  Table  7  gives  the  takeoff 
time  history  on  soil  and  Table  8,  the  takeoff  time  history  on 
a  paved  surface  thus  fulfilling  the  requirements  of  Problem  A  except 
for  the  percent  increase  in  takeoff  distance  which  is  merely  PI  *  total 
takeoff  distance  on  soil  minus  total  takeoff  distance  on  pavement 
times  100  divided  by  total  takeoff  distance  on  pavement  or  in  the 
actual  numbers  of  Problem  A: 

PI  -  3633-2240  X  100  (23) 

2240 


■  62% 

Table.  9  gives  the  results  for  Problem  B.  Note  that  for  any  speed,  the 
two  results  are  identical.  Finally,  Table  10  gives  results  for  Problem 
C.  A  drag  ratio  velocity  history  for  each  program  has  been  plotted  as 
Figure,  m  and  shows  the  afiv  t  of  simplification  on  the  firar  results. 
Notice  that  although  the  SOLDG  program  varies  from  the  results  of  S0LDG2 
and  TAKOFF,  it  varies  in  a  conservative  manner  from  a  design  point  of 
view. 

2.  Tables  of  Results  -  The  following  tables  are  computer  results  from 


Problems  A, 

B  and  C  presented 

in  Appendix  B. 
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Comparison  of  Sample  Problem  Results  Using  Landing 
Gear  Drag  Ratio  as  a  Function  of  Ground  Velocity 
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RITM  A  CONE  INOEX  OF  90.  ANO  CONSTANT  TAXI  VELOCITY 


TABLE  10  (Cont'd) 
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APPENDIX  D 


DATA  USED  TO  DESCRIBE  THE  BOEING  367-80  AIRCRAFT 


1.  Discussion  -  This  appendix  gives  the  detail  data  used  in  the  corre¬ 
lation  made  in  Section  IV.  Most  of  the  information  was  taken  from 
Reference  2.  A  number  of  items  were  not  specified,  however,  and  other 
sources  were  used  to  obtain  the  information.  The  takeoff  velocity 
required  for  lift-off  at  a  particular  gross  weight  was  obtained  from 
Reference  3,  and  is  presented  as  Table  11.  Another  problem  was  center 
of  gravity  variation  with  gross  weight.  It  was  convenient  in  the 
computer  program  to  use  the  center  of  gravity  as  number  of  feet  forward 
of  the  main  landing  gear  strut.  A  curve  of  Reference  2  was  converted 
to  the  convenient  form  of  Figure  1.1 2.  The  remaining  information  is 
presented  in  tabulated  form  below: 

367-80  Aircraft  Data 


Symbol  Used  in 

Item _ Computer  Program _ Value 


Number  of  Main  Landing  Gear  Tires 

•WMG 

8  per  strut 

Width  of  Main  Gear  Tire 

BM 

15  in. 

Diameter  ot  Main  Gear  Tire 

DM 

45  in. 

Main  Gear  Tire  Section  Height 

HTM 

10.75  in. 

Number  of  Nose  Landing  Gear  Tires 

WNG 

4 

Width  of  Nose  Gear  Tires 

BN 

10.65  in. 

Diameter  of  Nose  Gear  Tires 

DN 

30.7  in 

Nose  Tire  Section  Height 

HTN 

6.6  in 

Number  of  Main  Gear  Tandem  Wheel  Set* 

NTANDM 

4  per  strut 

Number  of  Nose  Gear  Tandem  Wheel  Sets 

NTANDN 

0 

Distance  Between  Nose  and  Main  Gear  Strut 

DWB 

44  Feet 

Distance  Between  Main  Cear  Strut  and  c.g. 

DCC 

See  Fig.  11 

Cross  Weight  of  Vehicle 

W 

120-200  KIP 

Wing  Area 

AW 

2821  Sq .  Ft 
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367-80  Aircraft  Data  (Cont'd) 


Item  Symbol  Used  in 

Item _ Computer  Program _ Value 


lift  Coefficient  During  Takeoff 

20*  Flap 

CL 

0.641 

Drag  Coefficient  During  Takeoff 

20*  Flap 

CD 

0.048 

Atmospheric  Density  Ratio 

SIG 

1 

Slope  of  Runway 

SLOPE 

0  rad 

Soil  Density 

RHOS 

2.6  lb  sec^ 

77“ 

Clay  Soil  Strength 

Cl 

80  to  200 

Sandy  Soil  Strength 

G 

20  to  150 

Paved  Surface  Friction  Coefficient 

RRC 

0.02 

Number  of  Engines 

NENGS 

4 

Thrust  o£  Engines 

XV EL,  YTH 

See  Figure 

Velocity  Required  for  Takeoff 

VT0 

See  Table  n 

Soil  Type 

NTYP 

Clay  and  Sand 

Load  Deflection  Curve  for  Main  Gear  Tire 
(46  x  16  Type  VII  26  PR  30  PS I 

Inflation) 

XFVM, 

WFM 

See  Figure  114 

Load  Deflection  Curve  for  Nose  Gear  Tire 
(32  x  11.50-15  12  PR  Type  VIII  at 

36  PSI  Inflation) 

XFVN , 

YDFN 

See  Figure  115 
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TABLE  11 

KC-135  TAKEOFF  SPEED 


SPEED-Knots 


Gross  Weight 

1000  Pounds 

20°  FlaD 

30°  Flap 

300 

184 

177 

290 

181 

174 

280 

177 

171 

270 

174 

163 

260 

171 

165 

250 

168 

162 

240 

164 

159 

230 

161 

155 

220 

157 

152 

210 

154 

148 

200 

150 

145 

190 

146 

141 

180 

143 

137 

170 

138 

134 

160 

134 

130 

150 

130 

125 

140 

126 

121 

130 

121 

117 

120 

116 

112 

110 

111 

107 

100 

100 

103 
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CG  Location  as  a  Function  of  Gross  Weight  for  Eoeing  367-80  Hi-FLota' 
Configuration 


213 


:ing  367-80  Engine  Takeoff  Thrust  as  a  Function 
True  Airspeed 


Figure  115  Boeing  367-80  Nose  Gear  Tire  Deflection 
Curve 
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APPENDIX  E 


DATA  USED  TO  DESCRIBE  THE  LOCKHEED  C-141  AIRCRAFT 


1.  Discussion  -  This  appendix  gives  basic  data  used  in  the  computer 
program  to  make  the  soil  takeoff  performance  assessment  of  Section  V. 
The  data  was  obtained  from  contractor  reports  and  discussion  with 
personnel  in  the  C-141  Systems  Program  Office  at  W-PAFB.  Some  of  the 
data  was  taken  directly  from  Reference  4  and  an  attempt  to  correlate 
takeoff  distance  on  paved  surfaces  was  made  with  the  results  contained 
therein.  Some  difficulty  was  encountered,  however,  since  thrust 
values  are  in  terms  of  Engine  Pressure  Ratio  (EPR)  and  not  as  installed 
thrust  as  required  by  the  present  program.  The  velocity  for  takeoff 
was  calculated  from  Equation  (24)  using  Lockheed  values  obtained 

from  the  stall  prevention  system  operation. 


V 


to 


(24) 


where:  V  =  Velocity  required  for  takeoff  (ft/sec) 

W  *  Vehicle  gross  weight  (lbs) 

S  =  Wing  Area  (sq  ft) 

P  *»  Atmospheric  density  (lb  sec^) 

ft* 

k  =  Takeoff  speed  safety  factor  (varies  from  1.1  to  1.5) 

CT  =  Lift  coefficient  for  stall  =  1.975 
ljMAX 

The  stall  speeds  were  increased  by  factors  of  1.1  to  1.5  from 
Figure  121  depending  on  the  thrust  to  weight  ratio.  The  results  are 
shown  in  figure  120.  Takeoff  criteria  used  in  this  study  was  to  obtain 
liftoff  speed.  *The  values  of  speed  and  gross  weight  used  in  the  analysis 
are  given  in  Table  12,  The  center  of  gravity  variation  with  gross  weight 
was  converted  to  distance  forward  of  the  main  landing  gear  strut  rs  in 
Appendix  D.  This  variation  is  shown  in  Figure  116.  Since  the  C-141 
aircraft  is  not  a  high  flotation  configuration  as  was  the  Boeing  367-80, 
much  higher  sinkage  and  drag  ratio  values  should  be  expected.  To  check  the 
effect  of  tire  pressure  on  the  performance  capability,  calculations 
were  made  for  tire  pressures  lower  than  the  rated  value  shown  in  the 
following  section.  These  tire  pressure  values  are  noted  on  specific 
curves  given  in  Section  V.  The  remaining  information  is  presented  in 
tabulated  form  on  c  following  page. 
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C-141  Aircraft  Data 


Symbol  Used  in 


_ Item _ Computer  Program 

Number  of  Main  Landing  Gear  Tires  WMG 

Width  of  Main  Gear  Tire  BM 

Diameter  of  Main  Gear  Tire  DM. 

Main  Gear  Tire  Section  Height  HTM 

Number  of  Nose  Landing  Gear  Tires  WNG 

Width  of  Nose  Gear  Tires  BN 

Diameter  of  Nose  Gear  Tires  DN 

Nose  Tire  Section  Height  HTN 

Number  of  Main  Gear  Tandem  Wheel  Sets  NTANDM 

Number  of  Nose  Gear  Tandem  Wheel  Sets  NTANDN 


Distance  between  Nose  and  Main  Gear  Strut  DWB 
Distance  Between  Main  Gear  Strut  and  c.g.  DCG 


Gross  Weight  of  Vehicle  W 

Wing  Area  AW 

Lift  Coefficient  During  Takeoff  25°  Flap  CL 

Drag  Coefficient  During  Takeoff  25°  Flap  CD 

Atmospheric  Density  Ratio  SIG 

Slope  of  Runway  SLOPE 

Soil  Density  RHOS 

Clay  Soil  Strength  Cl 

Sand  Soil  Strength  G 

Paved  Surface  Friction  Coefficient  RFC 

Number  of  Engines  NENGS 


_ Value _ 

4  per  strut 
15.68  in. 

42.85  in. 

12.22  in. 

2 

11.45  in. 

34.61  in. 

9.3  in. 

2  per  strut 
0 

53  ft. 

See  Figure  116 
140-316.1  KIP 
3228  Sq.  Ft. 
0.82 
0.055 
1 

0  rad 

2.6  ib  sec^/ft^ 
150  to  500 
45  to  200 
0.015 
4 
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C-141  Aircraft  Data  (Cont'd) 


Symbol  Used  In 


Item 

Computer  Program 

Value 

Thrust  of  Engines 

XVEL,  YTH 

See  Figure  117 

Velocity  Required  for 

Takeoff 

VTO 

See  Table  12 

Soil  Type 

NTYP 

Clay  and  Sand 

Load  Deflection  Curve 
(44x16  type  VII,  28 
Inflation) 

For  Main  Gear 
PR  130  PSI 

Tire 

XFVM, 

YDFM 

Set  Figure  118 

Load  Deflection  Curve 
(36x11  Type  VII,  22 

for  Nose  Gear 
PR,  140  PSI 

Tire 

XFVN, 

YDFN 

See  Figure  119 

Inflation) 
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TABLE  12 


C-141  TAKEOFF  SPEED 


Gross  Weight 

1000  lbs 

Speed  for 

Ft/Sec 

25°  Flap* 

Knots 

316.1 

238.28 

141.10 

300 

231.5 

137.09 

280 

223.68 

132.45 

260 

217.02 

128.51 

240 

209.9 

124.25 

220 

202.97 

120.19 

200 

195.49 

115.76 

180 

189.77 

112.37 

160 

189.87 

112.43 

140 

194.6 

115.25 

*  Speeds  correspond  to  handbook  lift-off  speed 
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PERCENT  MAC 


CC  LOCATION,  FT  FORWARD  OF  MAIN  GEAR 


TRUE  AIRSPEED  (KNOTS) 


7 


Figure  118  C— 141  Main  Gear  lire 
Deflection  Curve 


2J 


GROSS  WEIGHT,  1000  LB 


Takeoff  Speed  as  a  Function  of  Gross  Weight 
for  C-l4l  Aircraft 


THRUST/WEIGHT  AT  BRAKE  RELEASE 


Figure  121  Takeoff  Speed  as  a  Function  of  Thrust  to  Weight  Ratio 
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APPENDIX  F 


DATA  USED  TO  DESCRIBE  THE  LOCKHEED  C-5  AIRCRAFT 


1.  Discussion  -  This  appendix  gives  data  used  in  the  computer  program 
to  make  the  soil  takeoff  performance  calculations  which  appear  in  Section 
VI.  The  data  was  obtained  from  <1J  ;cussion  with  the  C-5  Systems  Program 
Office  and  from  contractor  reports.  The  stall  speed  was  determined  in 
the  usual  manner  by  equation  25  and  increased  by  a  factor  of  1.17  to 
obtain  takeoff  speed.  Table  13  shows  the  results  for  a  range  of  gross 
weights. 


V 


stall 


[  2W _ 

V  £  SCt 

V  lmax 


(25) 


'to 


=  K  V 


stall 


(26) 


where:  V  =  Velocity  required  for  takeoff  (ft/sec) 
to 

W  =  Vehicle  grots  weight  (lb) 

S  =  Wing  Area  (sq  ft) 

P  -  Atmospheric  density  (lb  sec^) 

k  =  Takeoff  speed  safety  factor  (1.17  all  weights) 

C^  =  Lift  coefficient  for  stall  =*  2.48  for  25°  flap 
MAX  iu  takeoff  configuration. 


The  center  of  gravity  variation  with  gross  weight  was  converted  to 
distance  forward  of  the  effective  main  gear  strut  and  is  given  as  Figure 
122.  This  concept  of  assuming  that  both  main  gear  trucks  are  mounted 
on  a  common  strut  was  chosen  for  convenience  since  the  program  was  designed 
assuming  only  ore  set  of  main  gear  struts.  This  approximation  gives 
higher  nose  gear  loads  and  equal  main  gear  loads.  The  overall  effect  of 
this  assumption  was  not  investigated  in  detail,  however,  it  was  felt  that 
results  would  be  representative  of  the  actual  performance  as  long  as  the 
correct  total  gross  weight  was  used.  The  remaining  data  is  presented 
in  tabular  form  on  the  following  page. 


C-5  Aircraft  Data 


Item _ 

Number  of  Main  Landing  Gear  Tlre8 

Width  of  Main  Gear  Tire 

Diameter  of  Main  Gear  Tire 

Main  Gear  Tire  Section  Height 

Number  of  Nose  Landing  Gear  Tires 

Width  of  Nose  Gear  Tires 

Diameter  of  Nose  Gear  Tires 

Nose  Tire  Section  Height 

Number  of  Main  Gear  Tandem  Wheel  Sets 

Number  of  Nose  Gear  Tandem  Wheel  Sets 

Distance  Between  Nose  and  Main  Gear 
Struts 

Distance  Between  Main  Gear  Strut 
and  c.g. 

Gross  Weight  of  Vehicle 
Wing  Area 

Lift  Coefficient  During  Takeoff,  25° 

Flap 

Drag  Coefficient  During  Takeoff,  25° 
Flap 

Atmospheric  Density  Ratio 
Slope  of  Runway 
Soil  Density 

Clay  Soil  Strength 


Symbol  Used  in 

Computer  Program _ Value 


WMG 

12  per  side 

BM 

17.16  in. 

DM 

47.94  in. 

HTM 

13.38  in. 

WNG 

4 

BN 

17.16  in. 

DN 

47.94  in. 

HTN 

13.38  in. 

NTANDM 

6  per  side 

NTANDN 

0 

DWB 

72.93  Ft. 

DCG 

See  Fig.  122 

W 

318.469-769  Kip 

AW 

6200  Ft2 

CL 

0.97 

CD 

0.0625 

SIG 

1 

SLOPE 

0  rad 

2 

RHOS 

2.6  lb  sec 
ft4 

Cl 

150  to  400 
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C-5  Aircraft  Data  (Cont*d) 


Item 

Symbol  Used  in 
Computer  Prosram 

Value 

Sand  Soil  Strength 

G 

45  to  200 

Paved  Surface  Friction  Coefficient 

RRC 

0.025 

Number  of  Engines 

NENGS 

4 

Thrust  of  Engines 

XVEL,  YTH 

See  Fig.  123 

Velocity  Required  for  Takeoff 

VTO 

See  Table  13 

Soil  Type 

NTYP 

Clay  and  Sand 

Load  Deflection  Curve  For 

Main  Gear  Tire  (49x17  Type  VII, 

26  PR,  70  to  130  PSI  Inflation) 

XFVM,  YDFM 

See  Figure  124 

Load  Deflection  Curve  for  Nose 

Gear  Tire  (Same  as  Main) 

XFVN,  YDFN 

See  Figure  124 

229 


TABLE  13 


C-5  TAKEOFF  SPEED 


Gross  Weight 
1000  Pounds 

800 

700 

600 

500 

400 

300 


Speed  -  Knots 
25°  Flap 

149.0 

139.0 

127.5 

114.5 

102.5 
89.0 
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Function  of  Gross  Weight  for  C-5  Aircraft 


SINGLE  ENGINE  THRUST  (1000  LB) 


0  20  40  60  80  100  120  140  160 


GROUND  VELOCITY  (KNOTS) 


Figure  123  C-5  Single  Engine  Thrust  me  a  Function 

of  Ground  Velocity 
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figure  124  Deflection  Curve  for  C-!>  Main  and  Noae 
Gear  Tires 
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